Guaranteeing Convergence of Iterative Skewed Voting Algorithms for Image Segmentation

In this paper we provide rigorous proof for the convergence of an iterative voting-based image segmentation algorithm called Active Masks. Active Masks (AM) was proposed to solve the challenging task of delineating punctate patterns of cells from fluorescence microscope images. Each iteration of AM consists of a linear convolution composed with a nonlinear thresholding; what makes this process special in our case is the presence of additive terms whose role is to "skew" the voting when prior information is available. In real-world implementation, the AM algorithm always converges to a fixed point. We study the behavior of AM rigorously and present a proof of this convergence. The key idea is to formulate AM as a generalized (parallel) majority cellular automaton, adapting proof techniques from discrete dynamical systems

Equiangular Tight Frames from Group Divisible Designs

An equiangular tight frame (ETF) is a type of optimal packing of lines in a real or complex Hilbert space. In the complex case, the existence of an ETF of a given size remains an open problem in many cases. In this paper, we observe that many of the known constructions of ETFs are of one of two types. We further provide a new method for combining a given ETF of one of these two types with an appropriate group divisible design (GDD) in order to produce a larger ETF of the same type. By applying this method to known families of ETFs and GDDs, we obtain several new infinite families of ETFs. The real instances of these ETFs correspond to several new infinite families of strongly regular graphs. Our approach was inspired by a seminal paper of Davis and Jedwab which both unified and generalized McFarland and Spence difference sets. We provide combinatorial analogs of their algebraic results, unifying Steiner ETFs with hyperoval ETFs and Tremain ETFs

Harmonic Equiangular Tight Frames Comprised of Regular Simplices

An equiangular tight frame (ETF) is a sequence of unit-norm vectors in a Euclidean space whose coherence achieves equality in the Welch bound, and thus yields an optimal packing in a projective space. A regular simplex is a simple type of ETF in which the number of vectors is one more than the dimension of the underlying space. More sophisticated examples include harmonic ETFs which equate to difference sets in finite abelian groups. Recently, it was shown that some harmonic ETFs are comprised of regular simplices. In this paper, we continue the investigation into these special harmonic ETFs. We begin by characterizing when the subspaces that are spanned by the ETF\u27s regular simplices form an equi-isoclinic tight fusion frame (EITFF), which is a type of optimal packing in a Grassmannian space. We shall see that every difference set that produces an EITFF in this way also yields a complex circulant conference matrix. Next, we consider a subclass of these difference sets that can be factored in terms of a smaller difference set and a relative difference set. It turns out that these relative difference sets lend themselves to a second, related and yet distinct, construction of complex circulant conference matrices. Finally, we provide explicit infinite families of ETFs to which this theory applies

The Bourgain-Tzafriri conjecture and concrete constructions of non-pavable projections

It is known that the Kadison-Singer Problem (KS) and the Paving Conjecture (PC) are equivalent to the Bourgain-Tzafriri Conjecture (BT). Also, it is known that (PC) fails for $2$-paving projections with constant diagonal $1/2$. But the proofs of this fact are existence proofs. We will use variations of the discrete Fourier Transform matrices to construct concrete examples of these projections and projections with constant diagonal $1/r$ which are not $r$-pavable in a very strong sense. In 1989, Bourgain and Tzafriri showed that the class of zero diagonal matrices with small entries (on the order of $\le 1/log^{1+\epsilon}n$, for an $n$-dimensional Hilbert space) are {\em pavable}. It has always been assumed that this result also holds for the BT-Conjecture - although no one formally checked it. We will show that this is not the case. We will show that if the BT-Conjecture is true for vectors with small coefficients (on the order of $\le C/\sqrt{n}$) then the BT-Conjecture is true and hence KS and PC are true

Harmonic Equiangular Tight Frames Comprised of Regular Simplices

An equiangular tight frame (ETF) is a sequence of unit-norm vectors in a Euclidean space whose coherence achieves equality in the Welch bound, and thus yields an optimal packing in a projective space. A regular simplex is a simple type of ETF in which the number of vectors is one more than the dimension of the underlying space. More sophisticated examples include harmonic ETFs which equate to difference sets in finite abelian groups. Recently, it was shown that some harmonic ETFs are comprised of regular simplices. In this paper, we continue the investigation into these special harmonic ETFs. We begin by characterizing when the subspaces that are spanned by the ETF\u27s regular simplices form an equi-isoclinic tight fusion frame (EITFF), which is a type of optimal packing in a Grassmannian space. We shall see that every difference set that produces an EITFF in this way also yields a complex circulant conference matrix. Next, we consider a subclass of these difference sets that can be factored in terms of a smaller difference set and a relative difference set. It turns out that these relative difference sets lend themselves to a second, related and yet distinct, construction of complex circulant conference matrices. Finally, we provide explicit infinite families of ETFs to which this theory applies

Harmonic Equiangular Tight Frames Comprised of Regular Simplices

An equiangular tight frame (ETF) is a sequence of unit-norm vectors in a Euclidean space whose coherence achieves equality in the Welch bound, and thus yields an optimal packing in a projective space. A regular simplex is a simple type of ETF in which the number of vectors is one more than the dimension of the underlying space. More sophisticated examples include harmonic ETFs which equate to difference sets in finite abelian groups. Recently, it was shown that some harmonic ETFs are comprised of regular simplices. In this paper, we continue the investigation into these special harmonic ETFs. We begin by characterizing when the subspaces that are spanned by the ETF\u27s regular simplices form an equi-isoclinic tight fusion frame (EITFF), which is a type of optimal packing in a Grassmannian space. We shall see that every difference set that produces an EITFF in this way also yields a complex circulant conference matrix. Next, we consider a subclass of these difference sets that can be factored in terms of a smaller difference set and a relative difference set. It turns out that these relative difference sets lend themselves to a second, related and yet distinct, construction of complex circulant conference matrices. Finally, we provide explicit infinite families of ETFs to which this theory applies

Auto-tuning Unit Norm Frames

Finite unit norm tight frames provide Parseval-like decompositions of vectors in terms of redundant components of equal weight. They are known to be exceptionally robust against additive noise and erasures, and as such, have great potential as encoding schemes. Unfortunately, up to this point, these frames have proven notoriously difficult to construct. Indeed, though the set of all unit norm tight frames, modulo rotations, is known to contain manifolds of nontrivial dimension, we have but a small finite number of known constructions of such frames. In this paper, we present a new iterative algorithm---gradient descent of the frame potential---for increasing the degree of tightness of any finite unit norm frame. The algorithm itself is trivial to implement, and it preserves certain group structures present in the initial frame. In the special case where the number of frame elements is relatively prime to the dimension of the underlying space, we show that this algorithm converges to a unit norm tight frame at a linear rate, provided the initial unit norm frame is already sufficiently close to being tight. By slightly modifying this approach, we get a similar, but weaker, result in the non-relatively-prime case, providing an explicit answer to the Paulsen problem: How close is a frame which is almost tight and almost unit norm to some unit norm tight frame

A Generalized Schur–Horn Theorem and Optimal Frame Completions

The Schur-Horn theorem is a classical result in matrix analysis which characterizes the existence of positive semi-definite matrices with a given diagonal and spectrum. In recent years, this theorem has been used to characterize the existence of finite frames whose elements have given lengths and whose frame operator has a given spectrum. We provide a new generalization of the Schur-Horn theorem which characterizes the spectra of all possible finite frame completions. That is, we characterize the spectra of the frame operators of the finite frames obtained by adding new vectors of given lengths to an existing frame. We then exploit this characterization to give a new and simple algorithm for computing the optimal such completion

Local Histograms and Image Occlusion Models

The local histogram transform of an image is a data cube that consists of the histograms of the pixel values that lie within a fixed neighborhood of any given pixel location. Such transforms are useful in image processing applications such as classification and segmentation, especially when dealing with textures that can be distinguished by the distributions of their pixel intensities and colors. We, in particular, use them to identify and delineate biological tissues found in histology images obtained via digital microscopy. In this paper, we introduce a mathematical formalism that rigorously justifies the use of local histograms for such purposes. We begin by discussing how local histograms can be computed as systems of convolutions. We then introduce probabilistic image models that can emulate textures one routinely encounters in histology images. These models are rooted in the concept of image occlusion. A simple model may, for example, generate textures by randomly speckling opaque blobs of one color on top of blobs of another. Under certain conditions, we show that, on average, the local histograms of such model-generated-textures are convex combinations of more basic distributions. We further provide several methods for creating models that meet these conditions; the textures generated by some of these models resemble those found in histology images. Taken together, these results suggest that histology textures can be analyzed by decomposing their local histograms into more basic components. We conclude with a proof-of-concept segmentation-and-classification algorithm based on these ideas, supported by numerical experimentation

Tremain Equiangular Tight Frames

Equiangular tight frames provide optimal packings of lines through the origin. We combine Steiner triple systems with Hadamard matrices to produce a new infinite family of equiangular tight frames. This in turn leads to new constructions of strongly regular graphs and distance-regular antipodal covers of the complete graph